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Abstract: Recent cosmological measurements favour additional relativistic energy density 
beyond the one provided by the three active neutrinos and photons of the Standard Model 
(SM). This is often referred to as "dark radiation", suggesting the need of new light states 
in the theory beyond those of the SM. In this paper, we study and numerically explore the 
alternative possibility that this increase comes from the decay of some new form of heavy 
matter into the SM neutrinos. We study the constraints on the decaying matter density 
and its lifetime, using data from the Wilkinson Microwave Anisotropy Probe, the South 
Pole Telescope, measurements of the Hubble constant at present time, the results from 
high-redshift Type-I supernovae and the information on the Baryon Acoustic Oscillation 
scale. We, moreover, include in our analysis the information on the presence of additional 
contributions to the expansion rate of the Universe at the time of Big Bang Nucleosynthesis. 
We compare the results obtained in this decaying matter scenario with those obtained with 
the standard analysis in terms of a constant A'cff . 
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1 Introduction 

In recent years an enormous progress has been made in the field of cosmology. Data on 
the Cosmic Microwave Background (CMB) has brought cosmology into a precision science 
and have revealed a Universe made by roughly 23% of Dark Matter (DM) and 72% of 
Dark Energy (DE). For a review on the physics of the CMB we refer to [1]. Furthermore, 
at present, recent analysis of cosmological data suggest a trend towards the existence of 
"dark radiation" (see for example [2-8] ) . Dark radiation alters the time of matter-radiation 
equality with a corresponding impact on the observed CMB anisotropics as well as an 
affect on the Large Scale Structure (LSS) distributions. The amount of dark radiation is 
usually parametrized using the parameter A'cff, that indicates the "effective number" of 
neutrino-like relativistic degrees of freedom. The value associated with the standard case 
of three active neutrino flavours was calculated in detail in Ref. [9] and was found to be 
iVSM ^ 3 046. 

The Wilkinson Microwave Anisotropy Probe (WMAP) collaboration found N^s = 
4.34+°!^ based on their 7-year data release and additional LSS data [2] at la in a ACDM 
cosmology. More recent measurements of the CMB anisotropy on smaller scales by the At- 
acama Cosmology Telescope (ACT) [10] and South Pole Telescope (SPT) [11] experiments 
seem to also favour a value of N^fi higher than predicted in SM. In a ACDM cosmology 



the current constraints on iVefj at 68% C.L at CMB time read [12]: 

AT^^MB ^ 4 34+0.86 WMAP7+BAO+i?o , (1-1) 

N^^'^ = 3.86 ± 0.42 WMAP7-FSPT+BAO-Fi?o , (1-2) 

^CMB ^ 3 gg ^ Q_4-^ WMAP7+ACT+SPT-FBAO+i?o , (1-3) 



where Hq refers to the constraint Hq = 74.2 it 3.6 km s found by the Hubble Space 
Telescope [13]. This evidence of dark radiation is robust under consideration of more 
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generalized cosmologies. For example in Ref. [3], N^g = A.Sb'^Q^^ was found in a global 
analysis including the data from cosmic microwave background experiments (in particu- 
lar from WMAP-7), the Hubble constant Hq measurement [13], the high-redshift Type-I 
supernovae[14] and the LSS results from the Sloan Digital Sky Survey (SDSS) data release 
7 (DR7) halo power spectrum [15]. The value refers to generalized cosmologies which de- 
part from ACDM models by allowing not only the presence of dark radiation but also dark 
energy with equation of state with u ^ —1, neutrino masses, and non- vanishing curvature. 

Independent information on the amount of radiation at earlier times is provided by its 
effect on the expansion rate of the Universe at the time of Nucleosynthesis. Faster expansion 
would lead to an earlier freeze-out of the neutron to proton ratio and would lead to a higher 
^He abundance generated during Big Bang Nucleosynthesis (BBN) [16]. Consequently the 
value of A^eff ^it the time of BBN can be constrained using primordial Nucleosynthesis yields 
of deuterium and helium. Old data and analysis reported a value of Agff ^.t BBN consistent 
with the standard model prediction: = 2.4 it 0.4 at 68% C.L. [17]. However recent 

results indicate, as well, that the relatively high "^He abundance can be interpreted in 
terms of additional radiation during the BBN epoch [18-20] (see [21] for a recent review). 
In particular in Ref. [20] it was found 

= 3.831°:^? (1.4) 

at 68% C.L. which suggests that AN^^^ < 1 at 95% C.L. at BBN time. This bound 
is reasonably independent of measurements on the baryon density from CMB anisotropy 
data and of the neutron lifetime input. 

The number of active neutrinos is constrained by measurements of the decay width of 
the Z boson [22] to be 2.984ib0.008. For this reason, several authors invoked the presence of 
a sterile neutrino to explain the above data [4, 23-26]. Furthermore, the possibility of one 
or two additional eV scale mass sterile neutrinos, brought to equilibrium with the SM ones 
via their mixing, could also account for some of the anomalies observed in short-baseline 
(SBL) neutrino experiments, that favour one or two sterile neutrinos [27-32]. 

Light sterile neutrinos account for dark radiation in the same amount at BBN and 
CMB times, i.e. in these scenarios = N^^^. Alternatively and in the light of the 

old BBN value for N^s [17], other models were proposed to explain an increasing on N^q 
at CMB time, while having Agfj compatible with three at BBN time. In Ref. [33], for 
example, the authors suggested that the increase in radiation could be explained by the 
decay of non-relativistic matter into relativistic states - see also [34-36] for a general study 
on dark radiation production from particle decays - while, in Ref. [37], it was proposed that 
non-thermal DM production, such as late-time decays of a long-lived state into DM and 
neutrinos and photons, could mimic an additional neutrino species. Recently, in Ref. [38], 
the authors discussed the possibility that a new particle of mass < 10 MeV, that remains 
in thermal equilibrium with neutrinos until it becomes non-relativistic, leads to a value of 
Agff that is greater than three. 

In this paper, we explore the possibility of generating the dark radiation without 
the need of additional light states besides the ones of the SM. In particular, we consider 
a scenario in which dark radiation consists of SM neutrinos which have been populated 
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in an extra amount by the decay of some heavy form of thermally produced matter as 
described in Sec. 2. One important constraint in our study is the requirement that the 
decay products, which constitute the dark radiation, are SM neutrinos. This is so because 
before neutrino decoupling any neutrinos produced by the decay would be brought to 
equilibrium by SM interactions and, in this way, the additional contribution from dark 
radiation to the expansion would be erased. 

We perform a global analysis of the relevant cosmological observables to determine the 
allowed range for the neutrino decaying matter density and its lifetime in this framework. 
Our results, presented in Sec. 3.2, show that this scenario can explain the tendency of the 
data to favour more radiation at the BBN and CMB times without the need of adding one 
or more sterile neutrinos nor any other new relativistic states to the three active ones of 
the SM. In particular we find that a lifetime of order Tdec ~ 10'^ s is favoured, implying a 
decay that happens during BBN time. We summarize our conclusions in Sec. 4. 

2 Decaying Matter 

Generically when discussing the possibility of decay matter as the origin of dark radiation 
one can think of two possibilities: in the first scenario the decaying matter is the main 
component of the DM itself, while in the second case the particle that decays is not the 
relic particle present nowadays in the Universe, but it is another particle species. 

The possibility of decaying DM (DDM) has been exhaustively explored in the litera- 
ture. CMB constraints on decaying warm DM have been derived in [39], specifically in the 
context of a Majoron DM. Cosmological constraints on the DM decay rate into neutrinos 
have been obtained in Ref. [40], using only Type la supernova data (SNIa) and the first 
CMB peak, see also [41-44] for more references. More recently, in Ref. [45], the authors 
extended the analysis using the full CMB anisotropy spectrum, large scale structure (LSS), 
Lyman-a data and weak lensing observations: the authors found a DM lifetime of the or- 
der r^g^ > 100 Gyr. Cosmological data, thus, tightly constrain the DM decay rate into 
neutrinos. This happens because DDM alters the time of matter-radiation equality, and 
the early integrated Sachs-Wolfe effect, with an increase in the first CMB peak. DDM 
also changes the late integrated Sachs- Wolfe effect, with a direct consequence on the CMB 
anisotropy spectrum at small multipoles, see [39] for a detailed explanation. Bounds on 
the DM lifetime can be obtained independently also from neutrino telescopes data. The 
Super-Kamiokande bounds on the DM flux from the Galactic Center, from the Earth and 
the Sun are given in [46] and, using these data, the constrain on the DM lifetime has been 
calculated in Refs. [47, 48]. 

Consequently since the DM decay rate into neutrinos is constrained to be much longer 
than the age of the Universe, the increase in the value of A'efr that we expect in the DDM 
scenario is, in general, tiny ^. 

Alternatively one can consider scenarios - denoted in the following as decM - in which 
the decaying state is different that the one dominantly present nowadays in the Universe as 

^Considering a DM lifetime that depends on time, as suggested in Ref. [49], it is possible to explain the 
excess in the radiation within the DDM scenario. 



-3- 



dark matter, such as those in Ref. [33-37]. In what fohows, we will consider that the decM 
is thermally produced and dominantly decays into neutrinos. We, moreover, restrict our 
study to the case of a non-relativisti decM at BBN time: its mass is greater than roughly 
10 MeV. 

2.1 Boltzmann Equations for Decaying Matter 

In the synchronous gauge, the line element is defined as: 

ds^ = a^{T){-dT^ + {6ij + hij)dx'dx^} , (2.1) 

where r is the conformal time, while t is the cosmological time {dt = air) dr), and a = R/Rq 
is the scale factor. The metric perturbation hij can be expanded in Fourier space as [50] 

hij{x, t) = J d^ke"^-^hkjh{k, r) + {kk - lSij)Gv{k, r)} , (2.2) 

with k = kk. We will use the notation defined above when writing the equations for the 
perturbations, i.e. the fields h{k,T) and r]{k,T). 

In the decM scenario, the set of Boltzman equations describing the background evolu- 
tion of the different components of the energy density of the Universe is complemented with 
two equations describing the evolution of the decM. and dark radiation density as [51, 52]: 

Pdec = -'iaHpdec — oTdec Pdec , (2-3) 
Pdr = -4:aHpdr + aVdecPdec, (2.4) 

where the over-dot denotes the derivative respect to conformal time r. The subscript "dr" 
indicates the dark radiation coming from the decM., while "dec" is the density of the decM. 
Tdec is the decM decay width. 

Correspondingly we derive the equations for the density fluctuations [50] generalized 
for the case of decM into SM neutrinos by adapting the equations for of DDM in [53-55] 
to the specifics of our case, in which the decaying products are SM neutrinos. We find: 



SnM = -^h, (2.5) 

= -^h- + aTdec^^ (Sdm - Sr) , (2.6) 

3 3 Pr 

Or = fc2 - - aTdec^eR , (2.7) 

ctr = \ - ^^^3 + + ^r)) - aT,J-^aR , (2.8) 

= [^^'-1 - + 1) - <^^dec—Fi , (2.9) 



where we have defined Pr = Pu + Pdn Pdm = Pc + Pdec, with piy the standard neutrino 
contribution, i.e. non coming from matter decay, and pc the density of standard cold dark 
matter. In Eq. (2.9) I > 3 and F2 = 2aR and we have used the conventions in Refs. [45, 50]. 
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Note that the density perturbation equation for the DM does not depend on the value of 
r^ec since both the background DM density and overdensity decay at the same rate [45] . 

We have implemented the previous equations for the background and the perturbations 
in the CLASS (Cosmic Linear Anisotropy Solving System) code [56]. We considered flat 
cosmology and an equation of state for the cosmological constant parameter given by 
a; = -L 

3 Numerical Analysis 
3.1 Cosmological Inputs 

In our analysis we include the results from the 7- year data of WMAP [2] on the temperature 
and polarization anisotropies, using the likelihood function as provided by the collabora- 
tion. The ACT [57] and SPT [11] experiments have probed higher multipole moments than 
WMAP. We implement the SPT data only in our analysis, since they give a bound in N^s 
of the same order as the one obtained considering both the ACT and the SPT data, see 
Eqs. (1.2), (1.3). We build the corresponding likelihood functions from the data, covari- 
ance matrix and window functions, introducing other three parameters in the analysis: the 
Sunyaev-Zel'dovich (SZ) amplitude, the amplitude of Poisson distributed point sources and 
the amplitude of clustered point sources. For these parameters, we used gaussian priors 
as given in Refs. [11, 58] and the templates of Ref. [11]. We fixed the foreground terms 
to be positive. In the following, with the notation "CMB" we will always refer to the 
combination of WMAP7 and SPT data. 

We introduce a Hubble parameter prior, based on the latest Hubble Space Telescope 
value [59]: Hq = 73.8 ± 2.4 km Mpc~^. This measurement of Hq is obtained from the 
magnitude-redshift relation of 240 low-z Type la supernovae at z < 0.1. 

We also include the luminosity measurements of high-z SNIa as given in Ref. [14]. This 
compilation, the "Constitution" set, consists of 397 supernovae and it is an extension of 
the previous sample, the "Union" set [60]. 

Finally we use the measurement of BAO scale obtained from the Two-Degree Field 
Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey Data Release 7 (SDSS 
DR7) [61]. We use as input data the two distance ratios at z = 0.2 and z = 0.35, with 
dz = T~ s{zd) / Dy [z) , where rgizd) is the comoving sound horizon at the baryon drag epoch 
and Dv{z) = [{l + z)'^D'^j^cz/H{z)Y/'^ ^ with Da the angular diameter distance and H{z) the 
Hubble parameter. We build the corresponding likelihood function using the covariance 
matrix as given in Ref. [61]. Since in the fitting procedure of [61], the value of dz is 
obtained by first assuming some fiducial cosmology {h = 0.72, Q^/i^ = 0.0223, ^rn = 0.25), 
we rescale the predictions for the BAO scale as explained in Ref. [62] . We do not introduce 
the information on the full power-spectrum of the SDSS DR7 survey [15], but we consider 
only the BAO measurement, since we will neglect the neutrino mass in our analysis. 

The analysis method we adopt is based on a Markov Chain Monte Carlo (MCMC) 
generator which employs the Metropolis-Hasting algorithm, see [3, 63] for more details. 
For convenience for the case of decM in the MC generator we use the parameters in the 
Table 1 for which we assume a flat prior. The parameters ng^r^As and F^gc/Gyr"^ are 
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respectively the scalar spectral index, the optical depth at reionization, the amplitude of 
scalar power spectrum at A; = 0.05 Mpc~^ and the decay rate into neutrinos. 

The first four parameters in Table 1 are related to the initial cosmological constant, 
baryon, DM and decM density as: 

p\ = nA-Hi (3.1) 

pi = • Hi • a"3 (3.2) 
pi = n,- Hi ■ (3.3) 

pj,,, = 0,ec-^o'-«"' (3.4) 

where H^ has a fixed value. The densities are all in units of [3c^/87rG], with c the speed 
of light and G the Newton's constant. Note that we have decided to use as MC 
parameter rather than Hq and that the parameter Hq is not equivalent to the Hubble 
parameter today. The "physical" Hubble parameter is calculated at each redshift using the 
Friedmann equation: 

H{z) = ^^^), (3.5) 

with i that runs over all the components of the Universe. After solving the Boltzmann 
equations, we obtain the present values {z = 0) of the Hubble constant Hq, of the baryon 
density and of the cold dark matter density Qc- More specifically, the latter parameter 
is defined as ilc = iPc{z = 0) + pdec{z = 0))/pcrit, but we will see that our results favour a 
negligible contribution of pdec to the dark matter density at present time. When presenting 
the results of our analysis, we will use these derived parameters (which we refer to as 
"analysis parameters") to better compare with the ACDM+A'eff case. 

At any time (or equivalently redshift z) the densities Pdec{z) and Pdr{z) can be trans- 
lated in a time dependent effective number of neutrinos depending on how this one is related 
with the observations. At BBN N^^ is determined by its contribution to the expansion 
rate of the Universe which is what affects the primordial abundances. Any form of energy 
density, relativistic or non-relativistic enters in this observable. Consequently, at BBN 

^ BBN ^ Pdec{zBT,^) +Pdr{zBBn) ^ BBN _ ^SM (3 5) 

Isi \Tl) PlV^BBMj 

with N'^^ = 3.046. To account for the constraints from BBN we impose in the analysis a 
prior on the value of Nf^^, Eq. (1.4) defined for definit eness, at the time at which BBN 
ends, ZBBN = rBBN/r° - 1 with Tbbn = 0.01 MeV (T^ =2.726 K). 

In CMB observables, the information on A'^eff arises mostly from its contribution to 
the determination of the the matter-radiation equality epoch. Thus, when translating our 
results in terms of N^q at CMB time, we will define 

AA^CMB = Pdr{zcMB) ^ ^3^^^ 

(I) (n) PlizCMB) 

with zcMB = 1100. 
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MC Parameters 



symbols 



Cosmological constant density 
Baryon density 
Dark Matter density 
Decaying Matter density 
Scalar spectral index 
Optical depth at reionization 

Amplitude of scalar power spectrum at A; = 0.05 Mpc 
Decay rate 



-1 



As 
rdec/Gyr 




r 



-1 



Table 1. Parameters used in the MC generation in the decM scenario. 



Furthermore, to implement thermal equilibrium condition on the dark radiation before 
neutrino decoupling, we impose 



with T^-dcc > 1 MeV . 

Using the above data and their theoretical predictions - obtained from the modified 
Boltzmann equations presented in the previous section - we construct the combined like- 
lihood function and the corresponding probability distribution function from which we 
obtain the one-dimensional and two-dimensional probability distributions, as described in 



For the sake of comparison, we also perform the corresponding analysis in the frame- 
work of a ACDM model with a constant AA'efr (but without imposing the condition 
Eq.(3.8)). In this case, the parameters used in the MC, for which flat priors are used, 
are the present Hubble constant Hq, baryon density r^b, cold dark matter density Vic and 
AA'cflf, together with ns,r, and As- 

3.2 Results 

The results of our cosmological fits are summarized in Tab. 2 and in Figs. 1 and 2. In 
Tab. 2, we report the best-fit values as well as the la and 95% C.L. range for the analysis 
parameters. We present the results for our model, the decM well as for the 

standard ACDM model with the addition of a fix N^s and with the implementation of 
the constraint on in both cases. We denote these two analysis as decM-l-BBN and 

ACDM-FTVeff+BBN. As seen from the table, in the decM-^BBN the set of data considered 
is better described with an amount of decaying matter with lifetime short enough to allow 
for a good fraction of it to have decayed into neutrinos at BBN time. 

In Fig. 1, we present the one-dimensional probabilities for the two analysis. For a 
better understanding, we report both the value of AA'cff at BBN and CMB time (which for 
ACDM-I-A'cff+BBN are identical by definition). Since the BBN data we are using favours 
values of AN^g < 1 at 95% C.L., adding the BBN information slightly changes the best-fit 
value for AiVcff in the ACDM-hiVcff+BBN analysis respect to the analysis of ACDM+Ncs 




= 3.046 



(3.8) 



Ref. [3]. 
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decM+BBN 


ACDM+iVcff+BBN 


Parameter 


best 


la 


95% 


best 


la 


95% 


Ho [km/s/Mpc] 


72.6 


+1.5 
-1.4 


+3.0 
-2.8 


73.0 


+ 1.5 
-1.5 


+2.8 

-3.0 


^bh?' X 100 


2.254 


+0.034 
-0.037 


+0.069 
-0.068 


2.258 


+0.032 
-0.037 


+0.065 
-0.070 




0.125 


+0.005 
-0.005 


+0.012 
-0.010 


0.127 


+0.006 
-0.006 


+0.011 
-0.012 


^Og{pdec/ P-t) at t = 10"^ s 


-4.61 


+0.61 
-0.73 


+0.92 
-1.7 


— 


— 


— 


Us 


0.973 


+0.009 
-0.009 


+0.018 
-0.018 


0.975 


+0.010 
-0.010 


+0.019 
-0.019 


T 


0.084 


+0.013 

-o!oi5 


+0.026 
-o!o26 


0.083 


+0.013 
-o!oi3 


+0.027 
-o!o24 


As X 10^ 


2.452 


+0.082 
-0.083 


+0.164 
-0.157 


2.449 


+0.075 
-0.083 


+0.155 
-0.159 


\og{Tdec/s) 


2.9 


+1.7 
-1.0 


+3.7 
-1.5 










0.50 


+0.30 
-0.19 


+0.58 
-0.42 


0.70 


+0.25 
-0.30 


+0.44 
-0.59 


AiVBBN 






< 0.90 


0.70 


+0.25 
-0.30 


+0.44 
-0.59 



Table 2. Best-fit values, 68% C.L. and 95% C.L. errors for the analysis parameters. 



(without the BBN information) and reduces the la range from AA'eff — 0.71i!^Q3g (without 
BBN) to AiVcff ^ 0.70t[]J^ (with BBN). 

We find that the preferred lifetime is Tdec ~ 10^ sec which means decaying slightly 
before or during BBN. This can be understood because at any time t the contribution of 
the decaying matter and the dark radiation produced in its decay to the energy density 
of the universe, parametrized as AA'^j|'^ and AA'^g respectively, can be approximated by 
integrating Eqs. (2.3) and (2.4): 



{t) 



AiV* 



t 

To 

Tde 

to 

11 

T 



Pdecjto) 
1 

2 Pde 



t 

^dec 



Erf 



t 

Tdec 



Tde 



to 



Pdecjto) 
2 p^(to) 



(3.9) 
(3.10) 
(3.11) 



where we have assumed that the decay occurs in a radiation dominated era, a oc \/t and 
we have chosen for normalization to = 10~^s (up to the factor Eq.(3.11) can be also 
obtained with the assumption of fast decay at t = Tdec [33]). 

At CMB time AA^^^^ is well approximated by Eq. (3.11) and the analysis constraints 
the product of Pdec{to) x ^Tdec so the density of decaying matter and its lifetime become 
strongly anticorrelated (see also Fig. 2). At BBN time 

'^i"" = AN^i^itBBN) + AAr,i^(tBBN) 

-) fliV'V^V'^^Erf 

7j \ aJ \to J p^ito) 2 [V Tde 



AN. 



tBBN 



(3.12) 
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Figure 1. Marginalized one-dimensional probability distributions for the analysis parameters. We 
present the results for the case of decM+BBN with solid lines. The dashed line shows the results 
for ACDM+A^cff+BBN. 

The present analysis favours N^^p^ not very different from AA^^^^^^ and Eq. (3.12) and 
Eq. (3.11) give comparable numerical results for r^ec ~ ^bbn- 

From the one-dimensional probability plots, it can also be inferred that the probability 
distributions for the parameter AA'gfr in the decM+BBN and in the model ACDM+A'^cfj+BBN 
are not equivalent. Besides the expected differences due to the fact that in the second sce- 
nario iVcfj is constant in time while in the first one it is not, there is a more subtle difference 
that is associated with the prior distribution assumed for the parameters in the two anal- 
ysis. In the decM model, A'eff is a derived parameter and it arises from a combination 
of two model parameters ildec and Fdec as given in Eq. (3.11) and Eq. (3.12). Let us in- 
troduce an auxiliary variable Z, function of ^dcc and T^cc^ which takes values along the 
A'cff=cnt direction. Using that the probability has to be invariant under reparametriza- 



-9- 




Figure 2. Two-dimensional probability distribution at the 68% and 95% C.L. for the decM pa- 
rameters {logipdec/Pi), log(Tdec/s)}, and {log(rdec/s), AN^fP^}. 



tions: T'{Qdcc,^dcc)d^decdTdec = T^iNcS, Z)dNesdZ. The flat prior for f^doc and Fdec that 
we are using in the analysis of the decM scenario means that P(Qdec; rdec)=cnt. Conse- 
quently, from Eqs.(3.9)-(3.11), we find that the jacobian for the change of variables leads 
to P{Nes,Z) oc 1/Neff. In other words, the effective prior in N^g, that corresponds to 
having a flat prior in {Odeo Tdcc}) pushes N^s to lower values. Furthermore, since Hq, 
ilfjh'^, VLch? and are correlated with the value of AiVcfj, this difference leads to the shift 
in the best-fit values for the parameters Hq, Qf^h'^, Qch'^ and in the two models seen in 
Tab. 2 and in Fig. 1. 

This also explains that the one-dimensional probability distribution for the parameter 
AA'^^'^ presents two peaks: one around AN^^ ~ and one around AA'^^^^ ~ 0.5. The 
first one is a consequence of the l/Ngg "initial" probability distribution in combination 
with the condition of vanishing dark radiation before neutrino decoupling, Eq. (3.8). The 
second is induced by the condition of having an amount of dark radiation at BBN, Eq. (1.4). 
Note, however, that the region below the first peak is much smaller than the region below 
the second peak. Indeed, the integral of fj'°^ dx Pi-dira{AN^^^) ~ 0.68, with Xmin = 0.28. 
This means that even if the height of the first peak is higher than the second one, a value 
of ANf^^ > 0.28 is favoured at roughly 68% C.L.. This is, in turn, consistent with the 
best-fit value of the lifetime Tdec, that is around BBN time (see right panel in Fig. 2). 

In the right panel of Fig. 2 we present the two-dimensional 68% and 95% C.L. credibility 
regions for the parameters log(rdec/s) and AA'^^^ for the decM-l-BBN analysis. In the 
figure, it is clearly visible that the regions are formed by two connected "islands" around 
the two favoured values of AA'^^^'^ . Also as expected the two parameters are anticorrelated: 
a shorter lifetime is associated with a value of AA^^^^^^ bigger than zero, while a lifetime 
bigger than 10^ will increase the value of N^^^ but not the value of A"^^^. 

The left panel of Fig. 2 shows the two-dimensional 68% and 95% C.L. credibility 
regions for the parameters {log{pdec/ P-y) and log(rrfec/s)}. Here we also see the expected 
anticorrelation between the total amount of decaying matter and its lifetime required to 
produce a certain amount of dark radiation (see Eq.(3.11)). The figure clearly shows how 
the best fit values of the analysis correspond to a decay during the BBN time while the 
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contribution of the favoured range of pdec implies that its contribution to the dark matter 
at present times is totally negligible. 

4 Conclusions 

In this paper, we have explored the possibility that dark radiation is not due to new rela- 
tivistic particles but it is formed by SM neutrinos whose density is increased by the decay 
of some heavy thermally produced particle. With this aim we have performed a global 
analysis of all the relevant cosmological data from the WMAP and SPT collaborations, 
measurements of the Hubble constant at present time, the results from high-redshift Type- 
I supernovae and the BAO scale. We have also included the information on additional 
radiation at the time of BBN as inferred from the relatively high ^He abundance deter- 
mination. We have concluded that the inclusion of the BBN information favours a decay 
that happens during BBN time with favoured lifetime T^ec between 10^ and 10^ seconds 
and leads to similar values of N^s at BBN and CMB times. We have discussed the dif- 
ference in the analysis between this scenario and the standard ACDM+A'cff associated 
with the different physical assumptions on the initial probability distribution of the model 
parameters. 

The discussion of a specific theoretical model in which this scenario is implemented is 
beyond the scope of this paper. Different possibilities have been suggested in the literature 
in the context of a decM that decays into radiation and that can be applied with some 
extents also to our case of decay into neutrinos. We refer to Refs. [33, 35] for different 
examples. 
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